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ABSTRACT 

The observed solar p-modc frequencies provide a powerful diagnostic of the in- 
ternal structure of the Sun and permit us to test in considerable detail the physics 
used in the theory of stellar structure. Amongst the most commonly used techniques 
for inverting such helioseismic data are two implementations of the optimally local- 
ized averages (OLA) method, namely the Subtractive Optimally Localized Averages 
(SOLA) and Multiplicative Optimally Localized Averages (MOLA). Both are con- 
trolled by a number of parameters, the proper choice of which is very important for 
a reliable inference of the solar internal structure. Here we make a detailed analysis 
of the influence of each parameter on the solution and indicate how to arrive at an 
optimal set of parameters for a given data set. 
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1 INTRODUCTION 

The observed solar p-mode oscillation frequencies depend 
on the structure of the solar interior and atmosphere. The 
goal of the inverse analysis is to make inferences about the 
solar structure given these frequencies. A substantial num- 
ber of inversions using a variety of techniques have been 
reported in the literature within the last decade (e.g. Gough 
& Kosovichev 1990; Dappen et al. 1991; Kosovichev 1993; 
Dziembowski et al. 1994; Basu et al. 1997). Two of the most 
commonly used inversion methods are implementations of 
the optimally localized averages (OLA) method, originally 
proposed by Backus & Gilbert (1968): the method of Multi- 
plicative Optimally Localized Averages (MOLA), following 
the suggestion of Backus & Gilbert, and the method of Sub- 
tractive Optimally Localized Averages (SOLA), introduced 
by Pijpers & Thompson (1992, 1994). Both methods de- 
pend on a number of parameters that must be chosen in 
order to make reliable inferences of the variation of the in- 
ternal structure along the solar radius. Most authors do not 
specify how these parameters are chosen or how a different 
choice would affect the solution. The goal of this work is to 
make a detailed analysis of the influence of each parameter 
on the solution, as a help towards arriving at an optimal set 
of parameters for a given data set. 

The adiabatic oscillation frequencies are determined 
solely by two functions of position: these may be chosen 
as density p and Ti = (d lnp/9 In p) a( j or as any other in- 
dependent pair of model variables related directly to these 
(e.g. Christensen-Dalsgaard & Berthomieu 1991). The solar 



p modes are acoustic waves that propagate in the solar in- 
terior and their frequencies are largely determined by the 
behaviour of sound speed c. Hence, it is natural to use c 
as one of the variables, combined with, e.g., p or Ti. The 
helium abundance Y is also commonly used, in combina- 
tion with p or p/p, p being pressure; this, however, requires 
the explicit use of the equation of state, incomplete knowl- 
edge of which could cause systematic errors (see Basu & 
Christensen-Dalsgaard 1997). In this work, we consider the 
inverse problem as defined in terms of sound speed and den- 
sity. 



2 LINEAR INVERSION TECHNIQUES 
2.1 The inverse problem 

Inversions for solar structure are based on linearizing the 
equations of stellar oscillations around a known reference 
model. The differences in, for example, sound speed c and 
density p between the structure of the Sun and the refer- 
ence model {5c? j<? ,Sp/p) are then related to the differences 
between the frequencies of the Sun and the model {8ui/ui) 
by 
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where r is the distance to the centre, which, for simplicity, 
we measure in units of the solar radius Rq . The index i num- 
bers the multiplets (n,l). The observational errors are given 
by e;, and are assumed to be independent and Gaussian- 
distributed with zero mean and variance of. The kernels 
K l 2 and K % 2 are known functions of the reference model. 
The term in F sur f (wj) is the contribution from the uncertain- 
ties in the near-surface region (e.g. Christensen-Dalsgaard & 
Berthomieu 1991); here Qi is the mode inertia, normalized 
by the inertia of a radial mode of the same frequency. 

For linear inversion methods, the solution at a given 
point ro is determined by a set of inversion coefficients a(ro), 
such that the inferred value of, say, 5c 2 /c 2 is 



Sc 2 



(2) 



From the corresponding linear combination of equations (|lj) 
it follows that the solution is characterized by the averaging 
kernel obtained as 



and also by the cross-term kernel: 
C(r ,r) = y>(r )^; c2 (r) , 



(3) 



(4) 



which measures the influence of the contribution from Sp/p 
on the inferred 5c 2 /c 2 . The standard deviation of the solu- 
tion is obtained as 



1/2 



(5) 



The goal of the analysis is then to suppress the contribu- 
tions from the cross term and the surface term in the linear 
combination in equation (j^), while limiting the error in the 
solution. If this can be achieved 



Sc 2 



IC(r , r)-4r(r)dr 
c 



(6) 



It is generally required that JC(ro,r) has unit integral with 
respect to r, so that the inferred value is a proper average 
of 5c 2 /c 2 : we apply this constraint here. Evidently, the res- 
olution of the inference is controlled by the extent in r of /C, 
the goal being to make it as narrow as possible. 

The surface term in equation ([!]) may be suppressed by 
assuming that F sur f can be expanded in terms of polynomials 

and constraining the inversion coefficients to satisfy 



CiMQrVAO*) =° > A = 0,1,..., A 



(7) 



(Dappen et al. 1991). As F aur f is assumed to be a slowly 
varying function of frequency, we use Legendre polynomials 
of low degree to define the basis functions ip\ . The maximum 
value of the polynomial degree, A, used in the expansion is 
a free parameter of the inversion procedures, which must be 
fixed. 

There are analogous expressions for the density inver- 
sion, expressing (Sp/p(ro)) in terms of the appropriate av- 
eraging kernel obtained as a linear combination of the mode 



kernels K % c2 , and involving a cross term giving the con- 
tribution from Sc 2 /c 2 . In the case of density inversion, an 
additional constraint is obtained by noting that the mass 
of the Sun is quite accurately known, the mass of the refer- 
ence model being usually fixed at this value; thus the density 
difference is generally constrained to satisfy 



4tt 



— (f) p(r) r 2 dr = . 
P 



(8) 



We have found that this constraint is important for stabi- 
lizing the solution. 

A number of different inversion techniques can be used 
for inverting the constraints given in equation (Q). We have 
used two versions of the technique of Optimally Localized 
Averages (OLA) (cf. Backus & Gilbert 1968) where the in- 
version coefficients are determined explicitly. 



2.2 SOLA Technique 

The aim of the Subtractive Optimally Localized Averages 
(SOLA) method (Pijpers & Thompson 1992, 1994) is to 
determine the inversion coefficients so that the averaging 
kernel is an approximation to a given target T(ro,r), by 
minimizing 

[JC(r ,r) - T{r ,r)] 2 dr + (3 f C 2 (r , r) f(r) dr 
o Jo 

+f 1& ~ 2 c 2 (r )(Ti , (9) 



subject to K. being unimodular. Here f(r) is a suitably in- 
creasing function of radius aimed at suppressing the sur- 
face structure in the cross-term kernel: we have used f(r) = 
(1 + r) 4 . Also, pi is a trade-off parameter, determining the 
balance between the demands of a good fit to the target 
and a small error in the solution; also, the quantity a 2 is the 
average variance, defined by 



M 



(10) 



M being the total number of modes. The second trade-off 
parameter j3 determines the balance between the demands 
of a well-localized averaging kernel and a small cross term. 
To suppress the influence of near-surface uncertainties, i.e., 
the term in F sur f, the coefficients are constrained to satisfy 
equation (Q). 

We have used target functions defined by 



T(ro,r) = Ar exp 



A(r 



r A(r ) 
2r 



(11) 



where A is a normalization constant to make the target 
unimodular. Thus the target function has its maximum at 
r — ro and has almost a Gaussian shape, except that it is 
forced to go to zero at r = 0. The target is characterized by a 
linear width in the radial direction: A(ro) = l\\c{ro) / c(rA), 
where ta is a reference radius; this variation of the width 
with sound speed reflects the ability of the modes to re- 
solve solar structure (e.g. Thompson 1993). We have taken 
ta = 0.2i?Q, and in the following characterize the width by 
the corresponding parameter Aa. 
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2.3 MOLA Technique 

In the case of Multiplicative Optimally Localized Averages 
(MOLA) method, the coefficients are found by minimizing 
1 pi 



K, (ro,r)J(ro,r)dr + (3 I C (ro,r) f(r)dr 



+/i<r 2 cf (r )q- 



2 

i : 



(12) 



where J(ro,r) is a weight function that is small near ro and 
large elsewhere: 



J(r ,r) = (r - r ) 2 



(13) 



This, together with the normalization constraint, forces K, 
to be large near ro and small elsewhere, as desired. As in 
equation (^) /(r) is included to suppress surface structure 
in the cross-term kernel. The quantity a 2 is defined by equa- 
tion (jlo|) . To suppress the influence of near-surface uncer- 
tainties, i.e., the term in F SUT f, the coefficients are again 
constrained to satisfy equation (Q). The MOLA technique 
is generally much more demanding on computational re- 
sources than is the SOLA technique because it involves anal- 
ysis of a kernel matrix which depends on the target ro; in 
the SOLA case, the corresponding matrix is independent of 
ro and hence need only to be analyzed once, for a given 
inversion case. 

2.4 Quality measures for the solution 

As seen from the previous sections, the inversions are char- 
acterized by the free parameters /j,, /3, A and Aa in the case 
of SOLA. These must be chosen to balance the relative im- 
portance of obtaining a well-localized average of the sound 
speed (or density) difference, minimizing the variance of the 
random error and reducing the sensitivity of the solution to 
the second function (i.e., the cross term) as well as to the 
surface uncertainties. 

The resolution of the inversion is characterized by the 
properties of the averaging kernel (eq. |^), which determine 
the degree to which a well-localized average of the underly- 
ing true solution can be obtained. Various measures of the 
width of K, have been considered in the literature. Here we 



measure resolution in terms of the distance A Q 



„(3) 



between the upper and lower quartile points of IC; these are 
defined such that one quarter of the area under /C lies to the 
left of rqu' and one quarter of the area lies to the right of 
rqu'. Furthermore, the location of r q u' and r$ , relative to 
the target location ro, provides a measure of any possible 
shift of the solution relative to the target. For the average of 
the solution to be well-localized, it is not enough that A qu 
be small: pronounced wings and other structure away from 
the target radius will produce nonlocal contributions to the 
average. As a measure of such effects in the SOLA case, we 
consider 



x(ro) 



[/C(r ,r)-T(r ,r)] 2 dr 



(14) 



which should be small (Pijpers & Thompson 1994). In the 
MOLA case, we introduce 



X(r ) 



K, (r ,r)dr+ / K, 2 (r ,r)dr 



(15) 
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Figure 1. Error correlation at ri = O.5i?0 for sound speed and 
density using SOLA and MOLA, based on the mode set described 
in Section 3; results are shown for a null surface term (F auI { = 0) 
and for A = 6, while the remaining inversion parameters have 
their default values (cf. Section 4)- 



where ta and rs are defined in such a way that the averag- 
ing kernel has its maximum at (ta + re)/2 and its FWHM 
is equal to (rs ~ ta)/^', again, a properly localized kernel 
requires that \' i s small. 

In a similar way, it is useful to define a measure C(ro) 
of the overall effect of the cross term: 



C(r ) = 



C 2 (r , r)dr 



(16) 



which should be small in order to reduce the sensitivity of 
the solution to the second function. 

It is evident that the overall magnitude of the error in 
the inferred solution should be constrained. However, Howe 
& Thompson (1996) pointed out that it is important to con- 
sider also the correlation between the errors in the solution 
at different target radii. This arises even if the errors in 
the original data are uncorrelated: the errors in the solu- 
tion at two positions are generally correlated, because they 
have been derived from the same set of data. The normalized 
correlation function which describes the correlation between 
the errors in the solution at r\ and at V2 is defined as: 



E(n,r 2 ) = 



■ 2 M)o~ 2 



1/2 



E 



c 2 (r 2 )* 2 ] 1/2 



(17) 



Howe & Thompson showed that correlated errors can intro- 
duce features into the solution on the scale of the order of 
the correlation-function width. 

Examples of correlation functions are shown in Fig. |l[ 
For sound-speed inversion, the error correlation generally 
has a peak at ri = r2 of width corresponding approximately 
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to the width of the averaging kernel (Fig. |l] top). For den- 
sity inversion, this peak at n = r-x is much broader than 
the averaging-kernel width. This is a consequence of the dif- 
ficulty in inferring density using acoustic-mode frequencies. 
There is also a region of strong anti-correlation (Fig. |l] bot- 
tom) . This is a result of applying the mass-conservation con- 
dition (eq. ^ since an excess of density in one part of the 
model has to be compensated by a deficiency in another. 



3 DATA AND MODELS 

The properties of the inversion depends on the mode selec- 
tion and errors in the data; the combination of mode selec- 
tion and errors is often described as the mode set, in contrast 
to the data set which in addition contains the data values. 
We have based the analysis on the combined LOWL + Bi- 
SON mode set described by Basu et al. (1997). Here the 
modes are in the frequency range 1.5-3.5 mHz, with degrees 
between and 99. This set in particular provides values for 
the standard errors o~i which to a large extent control the 
weights given to individual modes; here a 2 — 8.6 x 1CF 11 (cf. 
eq. id) . In some cases realizations of artificial data were con- 
sidered; these were obtained as differences between frequen- 
cies of the proxy and reference models, discussed below, with 
the addition of normally distributed random errors with the 
variances of the LOWL+BiSON mode set. 

Also, but to a far lesser extent, the inversion depends on 
the reference model. We have used Model S of Christensen- 
Dalsgaard et al. (1996) as our reference model. The model 
assumes that the Sun has an age of 4.6 Gyr. To construct 
artificial data for tests of the parameters for solar struc- 
ture inversion we adopted as a "proxy Sun" another model, 
of identical physical assumptions to those in Model S, but 
with the lower age of 4.52 Gyr. With no further modifica- 
tions, the "proxy Sun" would not include any surface uncer- 
tainties, and hence F 8ur f would have been zero. To provide a 
reasonably realistic model of this term, the "proxy Sun" in 
addition contained a near-surface modification based on a 
simple description of the effects of turbulent pressure on the 
frequencies, and calibrated to match the actual near-surface 
contribution to the difference between the solar frequencies 
and those of Model S (cf. Rosenthal 1998). 



4 THE CHOICE OF INVERSION 
PARAMETERS 

For a given mode set, the parameters controlling the inver- 
sion must be chosen in a way that, in an appropriate sense, 
optimizes the measures of quality introduced in Section 2.4. 
Needless to say, this places conflicting demands on the dif- 
ferent parameters, requiring appropriate trade-offs. Also, it 
is probably fair to say that no uniquely defined optimum 
solution exists. Here we have chosen what appears to be 
reasonable parameter sets, (cf. eqs |20| and ^l|) . The proce- 
dure leading to these choices is summarized in Section 4.4; 
however, we first justify them by investigating the effect on 
the properties of the inversion of modifications to the param- 
eters around these values. This is mostly done in terms of 
quantities such as error, error correlation, and kernel proper- 
ties which do not depend on the used data values; however, 
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Figure 2. On the left-hand side, the grey dots are QiSwi/uii after 
subtracting the fitted H\(u>i/L) and the black dots are the fits 
using A = 3 (top) and 6 (bottom). On right-hand side, residuals 
of the H2 fit are plotted. 



the effects are also illustrated by analyses of the artificial 
data defined in Section 3. 

The parameter A plays a somewhat special role, in that 
the suppression of the surface effects is common to both in- 
version methods (SOLA and MOLA) and to inversion for 5c 2 
and 8 p. For this reason we treat A separately, in Section 4.1. 
The response of the solution to the values of the remaining 
parameters depends somewhat on the choice of inversion 
method, and strongly on whether the inversion is for the 
sound speed or density difference. We consider sound-speed 
inversion in Section 4.2 and density inversion in Section 4.3. 



4.1 The choice of A 

Unlike the remaining inversion parameters the choice of the 
degree A used in the suppression of the surface term must 
directly reflect the properties of the data values; we base 
the analysis on the near-surface modification introduced in 
the artificial data according to the procedure of Rosenthal 
(1998) (very similar results are obtained for solar data). 
To determine the most appropriate value of A we consider 
the frequency-dependent part of the frequency differences. 
This is isolated by noting that according to the asymptotic 
theory the frequency differences satisfy (e.g. Christensen- 
Dalsgaard, Gough & Thompson 1989) 



5u>i 



Hi 



(18) 



with L — 1 + 1/2, where / is the degree of mode i. Here Si is a 
scaling factor which in the asymptotic limit is proportional 
to Qi and the slowly varying component of Hai^i) corre- 
sponds to the function F sur f in the asymptotic limit. Thus, 
by fitting a linear combination of Legendre polynomials to 
H 2 : 
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4.2 Parameters for sound-speed inversion 




As reference we use what is subsequently determined to be 
the best choice of parameters: 



Figure 3. Trade-off diagram for sound-speed inversion at two 
different radii using SOLA (left) and MOLA (right). Values of 
/i and Aa increase from top to bottom; and f3 from bottom to 
top. In the SOLA case, /i varies from 10~ 5 to 10" 1 , (3 from 0.1 
to 1000 and Aa from 0.04 to 0.08. In the MOLA, ^ varies from 
10 — 6 to 10 — 2 and /3 from 0.1 to 1000. Parameters not explicitly 
mentioned have their reference values (cf. eq. [g^ ). Note that /i = 
10 _1 in the SOLA case leads to a large width when ro = 0.06-Rq 
and is outside the plot. 



H'2(cui) ~ a\P\(tUi) 



(19) 



we can determine the appropriate value of A for any given 
data set. In practice, we make a non-linear least-squares fit 
to a sum of two linear combinations of Legendre polynomi- 
als, in to/L and lu to SiSuji/uji, using a high A (A = 16). 
Then we remove Hi from SiSui/ui and fit now a single lin- 
ear combination of Legendre polynomials in u, looking for 
the smallest value of A that provides a good fit (Fig. |^). On 
this basis, we infer that A = 6 provides an adequate repre- 
sentation of the surface term; we use this as our reference 
value in the following. The solar data considered by Basu 
et al. (1997) have a similar behaviour, and A = 6 is also an 
appropriate choice in that case. 

The constraints imposed by equation (0) do not depend 
explicitly on the target location; hence it is reasonable that 
they introduce a contribution to the errors in the solution 
that varies little with ro, leading to an increase in the error 
correlation. This is confirmed in the case of sound-speed in- 
version by the results shown in the top panel of Fig. [I]. (We 
note that, in contrast, for density inversion the correlation 
decreases somewhat with increasing A; we have no expla- 
nation for this curious behaviour, but note that the density 
correlation is in any case substantial.) 
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(20) 



Effects on the quality measures of varying the parameters 
around these values are illustrated in Figs |§| [l|, ^| and ^| 
in addition, Fig. ^ shows results of the analysis of artifi- 
cial data (cf. Section 3), and Fig. |§] illustrates properties 
of selected averaging kernels. Throughout, parameters not 
explicitly mentioned have their reference values. 



4-2.1 Choice of fi 

Generally known as the trade-off parameter, /i must be de- 
termined to ensure a trade-off between the solution error 
(eq. ^) and resolution of the averaging kernel. This is typi- 
cally illustrated in trade-off diagrams such as Fig. ^, showing 
the solution error against resolution (here defined by the sep- 
aration A qu between the quartile points) as p varies (circles) . 
As n is reduced, the solution error increases; the resolution 
width generally decreases towards a limiting value which, in 
the SOLA case, is typically determined by the target width 
A(ro). On the other hand, for larger values of /x, there is 
a strong increase in the width, with a corresponding very 
small reduction in the solution error. 

The behaviour in the trade-off diagram depends on the 
target radius ro considered, the risk of a misleading solution 
being particularly serious in the core or near the surface if 
H is too large. Thus it is important to look at the trade-off 
diagram at different target radii. This is illustrated in Fig. ^, 
where solution error (lower panels) and the location of the 
quartile points relative to the target radius (upper panels) 
are plotted. It is evident that the error increases markedly 
towards the centre and surface, particularly in the SOLA 
case. In addition, the averaging kernels get relatively broad 
and there is a tendency that they are shifted relative to the 
target location, particularly near the centre. 

Note that the results plotted in Fig. ^ use values of 
[i such that the solution error given by SOLA and MOLA 
techniques are similar. Even for small values of [i, the MOLA 
averaging kernels do not penetrate as deep into the core as do 
the SOLA averaging kernels. The resolution of the averaging 
kernel is more sensitive to the choice of [i using MOLA than 
using SOLA (cf. Fig. §). 

In addition to the error and resolution, we also need to 
consider other properties of the solution. 'Global' properties 
of the averaging kernels, measured by \ or x' are illustrated 
in Fig. [B], together with the integrated measure C of the 
cross-talk. For larger values of n, these quantities increase, 
particularly near the surface. The strong increase in x' (and 
X) at large target radii is due to the presence of a depres- 
sion in the averaging kernel near the surface that increases 
quickly with ro (cf. Fig. ^), especially for MOLA. The influ- 
ence of /j, on the error correlation is illustrated in Fig. ^; the 
correlation evidently increases with decreasing fi, together 
with the solution error. Note that the error correlation in- 
creases with the target radius for any choice of parame- 
ters (see also Rabello-Soares, Basu & Christensen-Dalsgaard 
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Figure 4. Resolution of averaging kernels for sound-speed inversion using SOLA (left) and MOLA (right), illustrated by the locations 
r qu ~ r and riu ~ r of the upper and lower quartile points relative to the target radius (top), and solution error (bottom), as a function 
of target radius, for different values of fi. 
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Figure 5. Variation of C and \ (x! * n the MOLA case), for sound-speed inversion. In the SOLA case, fj, = 1CT 5 and 10" 3 , (3 = 0.1 
and 10, and Aa = 0.04 and 0.08; in the MOLA case, fi = 10~ 6 and 10 -4 and f} = 0.1 and 10. The continuous thick lines use the 
reference choice of parameters ( cf. eq. [sfjj; parameters not explicitly mentioned have their reference values. In the SOLA case, a larger 
value of Aa has no visible effect on \, while C is insensitive to the change in Aa- 
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Figure 6. Variation of error correlation with target radius using SOLA (left) and MOLA (right) for sound-speed inversion. In the 
SOLA case, fx = 10 -5 and 10 -3 , /3 = 0.1 and 10, and A A = 0.05 andj).07. In the MOLA case, fx = 10 -6 and 10 -4 and (3 = 0.1 and 
10. Again, the continuous thick lines use the 'optimal' values (cf. eq. \20\ ), and parameters not explicitly mentioned have their reference 
values. The same symbols are used as in Fig. H. 
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Figure 7. Solution (5c 2 /c 2 ) versus radius for different /i using SOLA (left) and MOLA (right), for the artificial data described in 
Section 3 and a single realization of errors. Vertical bars show one-sigma errors in the inferred solution, and are plotted at the second 
quartile point r^ , while the horizontal bars extend between the first and third quartile points r^ and r^l . Here r^ , r^} and r^ are 
defined such that 25%, 50% and 75% of the area under the averaging kernel lies to the left of them, respectively. The middle row uses 
the 'optimal' choice of parameters. The thick line is the difference between the theoretical models used. The thin line is the solution of 
the inversion without adding errors to the calculated eigenfrequencies. 
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Figure 8. Top: Averaging kernels for SOLA sound-speed inversion at different target radii ro and different Aa- Note that for a given 
Aa the width of the averaging kernel is reduced as the target radius increases. This is because the width changes with the sound speed 
(Section 2.2). Bottom: Detail of the averaging kernels showing their wings for two target radii and two values of Aa- Near the surface, 
there is a negative lobe resulting from the fact that our mode set has modes only for I < 99. Closer inspection of the computed kernels 
indicates that we do not have sufficient information to go further than ro ~ 0.95. 



1998). It is slightly smaller using SOLA than MOLA when 
their solution errors are similar. 

Finally, Fig. ^ shows the inferred solutions for the 
sound-speed difference obtained from the artificial data de- 
scribed in Section 3, for three values of n, and compared 
both with the true 8c? j<? and the solution inferred for error- 
free data. The behaviour of the solution generally reflects the 
properties discussed so far. In particular, it should be noticed 
that the solution for the data with errors is shifted system- 
atically relative to the solution based on the error-free data, 
reflecting the error correlation, most clearly visible in the 
outer parts of the model; this behaviour illustrates the care 
required in interpreting even large-scale features in the solu- 
tion at a level comparable with the inferred errors. We also 
note that even the inversion based on error-free data shows 
a systematic departure from the true solution, particularly 
near the surface. This appears to be a residual consequence 
of the imposed near-surface error, exacerbated by the lack 
of high-degree modes which might have constrained the so- 
lution in this region. We have checked this by considering 
in addition artificial data without the imposed near-surface 
error (cf. Section 3). 



4.2.2 Choice of (3 

The importance of f3 is seen most clearly in Fig. [s| in terms 
of the properties of the averaging kernels and the cross term: 
as desired, increasing j3 reduces the importance of the cross 
term as measured by C, but at the expense of poorer aver- 
aging kernels, as reflected in \ an d x' ■ It should be noticed, 
however, that the choice of mainly affects the solution in 
the core and near the surface, while it has little effect in the 
intermediate parts of the solar interior. We also note that the 
error in the solution increases with increasing (3 (cf. Fig. ^), 
as does the error correlation (see Fig. ^). Thus the choice 
of /3 is determined by the demand that the cross term be 
sufficiently strongly suppressed, without compromising the 
properties of the averaging kernels and errors. 

4.2.3 Choice of A A 

The SOLA technique has an additional parameter: the width 
of the target function at a reference radius (see eq. |ri| ). The 
aim of SOLA is to construct a well-localized averaging kernel 
that will provide as good a resolution as possible. As illus- 
trated in Fig. Aa ensures a trade-off between averaging- 
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0.05 0.10 0.05 0.10 

Figure 9. Trade-off diagram for density inversion at two differ- 
ent radii using SOLA (left) and MOLA (right). Values of /i and 
Aa increase from top to bottom: /x varies from 10 -7 to 10 — 3 
(SOLA) and from 10~ 9 to 10~ 5 (MOLA), whereas A a varies 
from 0.04 to 0.08. Parameters not explicitly indicated correspond 
to the best choice (cf. eq. [j"^). Note that fi = 10 — 3 in the SOLA 
case corresponds to a large width, outside the plot. 

kernel resolution (taking into account also the deviation x 
from the target) and the solution error. 

The effect of Aa on the averaging kernels is illustrated 
in Fig. ^. Evidently, for high Aa the solution is smoothed 
more strongly than at low Aa- However, if Aa is too small, 
the averaging kernel starts to oscillate; even more problem- 
atic is the presence of an extended tail away from the target 
radius since it introduces a non-zero contribution from radii 
far removed from the target. As in the case of /x, the error 
increases with increasing resolution when Aa is reduced. On 
the other hand, the error correlation decreases with decreas- 
ing Aa due to the stronger localization of the solution (cf. 
Fig. ^ and develops a tendency to oscillate. 

We finally note that C is almost insensitive to Aa. 

4.3 Density inversion 

As reference we use what is subsequently determined to be 
the best choice of parameters: 

SOLA : A = 6 , fi = 1(T 5 , f3 = 10 , A A = 0.06 ; 
MOLA : A = 6, /x = 10~ 7 , P = 50 . (21) 

Effects on the quality measures of varying the parameters 
around these values are illustrated in Figs |{| [lo| [ll] and |l2| 
in addition, Fig. |l3| shows results of the analysis of artificial 
data (cf. Section 3), and Figs |l4| and [l5| illustrate properties 
of selected averaging kernels. Throughout, parameters not 
explicitly mentioned have their reference values. 



In the case of density inversion, we have found that the 
cross term is generally small, with little effect on the so- 
lution. In addition, the sound-speed difference between the 
Sun and calibrated solar models is typically small in the 
convection zone (e.g. Christensen-Dalsgaard & Berthomieu 
1991), further reducing the effect of the sound-speed contri- 
bution in the density inversion. As a result, the effect of the 
value of (3 on the properties of the inversion is very modest, 
although a value of /3 in excess of 1 is required to suppress 
the remaining effect of the cross term. Hence, in the follow- 
ing we do not consider the effect of changes to /3, or the 
behaviour of C. 



4-3.1 Choice of /x 

As for sound-speed inversion, the trade-off parameter fi must 
be determined to ensure a trade-off between the solution er- 
ror and resolution of the averaging kernel (Fig. ^| - circles). 
Its behaviour is very similar to that in the sound-speed case 
(Fig. 0). As /x is reduced, the solution error increases but 
the resolution width cannot get smaller than a certain value, 
which, in the case of SOLA, is the target width A(ro). On 
the other hand, for larger values of /x, there is a strong in- 
crease in the width, with a corresponding very small reduc- 
tion in the solution error. 

The dependence of the trade-off on target radius ro is 
illustrated in Fig. If /x is too large, one may get a mislead- 
ing solution especially in the core. As for sound-speed inver- 
sion, the averaging-kernel resolution using MOLA is more 
sensitive to variations in /j, than using SOLA. 

For larger values of /x, beside the increase in the 
averaging- kernel width there is an increase in \ an d x' (°f- 
Fig. hll). The bump in \ an d X around ro — 0.3 for a large 
/x is due to a "shoulder" in the averaging kernel that ap- 
pears at these target radii, as illustrated in the case of SOLA 
in Fig. [l4| The error correlation, illustrated in Fig. [l^, in- 
creases somewhat with increasing /x; as already noted, the 
error correlation changes sign and is of much larger magni- 
tude for density than for sound speed, probably as a result 
of the mass constraint (eq. ^). 

Examples of inferred solutions, for the artificial data 
defined in Section 3, are shown in Fig. jl3|, together with 
the true model difference and the difference inferred from 
error-free frequency differences. For a large value of /x the 
averaging kernel is not well localized and hence affects the 
solution (bottom panel); this is true also for the solution 
based on error-free data. For smaller /x the error clearly in- 
creases; also, particularly in the SOLA case, the effect of the 
error correlation is evident. Thus the solution again deteri- 
orates. This illustrates how the correlated errors can intro- 
duce features into the solution, showing the importance of 
limiting the error correlation. 

4.3.2 Choice of A A 

As for the sound-speed inversion, Aa ensures a trade-off 
between averaging-kernel resolution and the solution error 
in the SOLA technique (cf. Fig. |^). As before, if we choose 
Aa too small, the averaging kernel is poorly localized and 
it starts to oscillate; this is reflected in the departure of x 
from the target (see Fig. [H]) and illustrated in more detail, 
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Figure 10. Resolution of averaging kernels for density inversion using SOLA (left) and MOLA (right), illustrated by the location of 
the upper and lower quartile points relative to the target radius (top), and solution error (bottom), as a function of target radius, for 
different values of p.. 
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Figure 11. Variation of x (x' * n the MOLA case) for density inversion. In the SOLA case, p, = 10 — 6 and 10 — 4 , and Aa = 0.04; in 
the MOLA case, fi = 10 -8 and 10 -6 . The continuous thick line uses the best choice of the parameters (cf. eq. ^\). All the inversions 
use the best choice of parameters, except when explicitly indicated. A larger Aa value cannot be distinguished from the best choice. 



for selected target radii, in Fig. ttq. On the other hand, for 
large Aa, the solution is smoothed relative to the one for 
small A a (due to the low resolution). The effect on the error 
correlation of changes in A a is very modest, although for 
small Aa there is a tendency for oscillations (cf. Fig. [l2f), as 
was also seen for sound-speed inversion. 



4.4 Summary of the procedure 

As a convenience to the reader, we briefly summarize the se- 
quence of steps which we have found to provide a reasonable 
determination of the trade-off parameters for the SOLA and 



MOLA methods for inversion for the corrections Sc 2 and 8p 
to squared sound speed and density: 

• The parameter A is common to both inversion methods 
and to inversion for Sc 2 and 5p. Unlike the other parameters, 
it must directly reflect the properties of the data values, in 
terms of the ability to represent the surface term (cf. Section 
4.1). It is also largely independent of the choice of the other 
parameters: fi, (5 and Aa; thus its determination is a natural 
first step. It should be noted, however, that the choice of A 
has some effect on the error correlation (cf. Fig. [l]) , generally 
requiring that A be kept as small as possible. 

• The second step is the determination of fb, whose value 
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case, ft = 1CT 6 and 1CT 4 , and A A = 0.04. In the MOLA case, fi = 10 8 and 10 6 . Again, the continuous thick line uses the 'optimal' 
values (cf. eq. \2l\ ). All the inversions use the best choice of parameters, except when explicitly indicated. A larger Aa value cannot be 
distinguished from the best choice. The symbols are as in Fig. \lA 



SOLA 



MOLA 



o 

X 





-50 
-100 



-50 
-100 



-50 
-100 




Figure 13. Solution (Sp/p) versus radius for different fi using SOLA (left) and MOLA (right). The middle row uses the best choice of 
parameters. The thick line is the difference between the theoretical models used. The thin line is the solution of the inversion without 
noise to the calculated eigenfrequencies. Note how the error correlation is responsible for a bad solution when using low fi. 
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Figure 14. Averaging kernels for density inversion at different 
target radii ro and different fj, using the SOLA method. The aver- 
aging kernels for ro = 0.3-Rq and OARq are displaced vertically 
for clarity. The thick lines use the best choice of the parameters 
(eg. [jjj. The full and empty circles use p, = 10 -4 and 10~ 6 , 
respectively. 



is the most critical to achieve a good solution. As described 
in Sections 4.2.1 and 4.3.1, it must be determined to ensure 
a trade-off between the solution error and resolution of the 
averaging kernel (Figs ^| and [] - circles) at representative 
target radii ro (see also Figs ^ and |l^). In addition, we need 
to consider the broader properties of the averaging kernels 
as characterized by \ (SOLA) or x[ (MO LA) and the cross- 
talk quantified by C (Figs ft and hit), as well as the error 
correlation (Figs H and ha). 

• The next step is to find /3 which is determined by the de- 
mand that the cross term be sufficiently strongly suppressed, 
without compromising the properties of the averaging ker- 
nels and errors. Its effect on the properties of density inver- 
sion is very modest. 

• Finally, the SOLA technique has an additional parame- 
ter: Aa, which ensures a trade-off between averaging- kernel 
resolution and the solution error. Aa is typically decreased 
until the averaging kernels are poorly localized and start to 
oscillate (Figs [i] and |l5|) which is reflected in the departure 
of x from the target. 

After this first determination of /x, /3 and possibly Aa, 
we go back to step number 2, determining /1 using now the 
new values of (3 and Aa- The procedure obviously requires 
initial values of (3 and (for SOLA) Aa: we suggest /3 = 10 
and A a = 0.06 or larger. 



Although the measures of quality of the inversion are 
essentially determined by the mode set (modes and errors), 
the determination of the parameters must also be such as to 
keep the solution error sufficiently small to see the variations 
in the relative sound-speed or density differences, which in 
the case of solar data and a suitable reference model could 
be as small as 10 -3 and 5 x 10 -2 respectively. Furthermore, 
to enable comparison of the solution of the inversion of two 
different data sets, the solution errors should be similar. To 
obtain an impression of the quality of the solution and the 
significance of inferred features, we also strongly recommend 
analysis of artificial data for suitable test models, including 
comparison of the inferred solutions for the selected param- 
eters with the true difference between the models and with 
the solution inferred for error-free data (see Figs and |l^) . 



5 CONCLUSION 

Appropriate choice of the parameters controlling inverse 
analyses of solar oscillation frequencies is required if reli- 
able inferences are to be made of the structure of the solar 
interior. This choice must be based on the properties of the 
solution, as measured by the variance and correlation of the 
errors, by the resolution of the averaging kernels and by the 
influence of the cross-talk. We have considered a mode set 
representative of current inverse analyses and investigated 
the properties of the inversion, as well as the solution cor- 
responding to a specific set of artificial data. By varying 
the parameters we have obtained what we regard as a rea- 
sonable choice of parameters (cf. eqs ^ and ^l|); this was 
verified by considering in some detail the sensitivity of the 
relevant measures of the quality of the inversion to changes 
in the parameters. The analysis also illustrated that an un- 
fortunate choice of parameters may result in a misleading 
inference of the solar sound speed or density; furthermore it 
became evident that the correlation between the error in the 
solution at different target location plays an important role, 
even for our optimal choice of parameters, and hence must 
be taken into account in the interpretation of the results (see 
also Howe & Thompson 1996). 

The meaning of the parameters is evidently closely re- 
lated to the precise formulation of the inverse problem. For 
example, it would be possible to introduce weight functions 
in the integrals in equations (^) and (|l2|), to give greater 
weight to specific aspects of the solution. The need for such 
refinements is suggested by the fact that the properties of 
the solution depends rather sensitively on the target loca- 
tion ro. More generally, it is likely that the best choice of 
parameters may depend on the target location, further com- 
plicating the analysis and (particularly in the SOLA case) 
increasing the computational expense. 

The procedure adopted here is evidently somewhat ad 
hoc, although we have attempted a logical sequence in the 
order in which the parameters were chosen. A more sys- 
tematic approach, making use of objective criteria, would 
in principle be desirable. However, even in the consider- 
able simpler case of inversion for a spherically symmetric 
rotation profile, characterized essentially just by the param- 
eters /x and possibly Aa, such objective determination of 
the parameters has so far met with little success in prac- 
tice (see, however, Stepanov & Christensen-Dalsgaard 1996 
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Figure 15. Top: Averaging kernels for SOLA density inversion at different target radii ro and different Aa- Note that for a given Aa 
the width of the averaging kernel is reduced as the target radius increases. This is because the width changes with the sound speed (see 
Section 2.2). Bottom: Detail of the averaging kernels showing their wings for two target radii and two values of Aa- As in the case of 
the sound-speed inversion, the restriction of our mode set to I < 99 causes a negative lobe near the surface. Closer inspection indicates 
that we do not have sufficient information to go further than ro ~ 0.92. 



and Hansen 1996). On the other hand, it is far from obvious 
that an objectively optimal solution to the inverse problem 
exists, for a given data set: the best choice of parameters 
may well depend on the specific aspects of the solar inte- 
rior that are being investigated. It is important, however, 
that the error and resolution properties of the solution be 
kept in mind in the interpretation of the results; indeed, the 
immediate availability of measures of these properties is a 
major advantage of linear inversion techniques such as those 
discussed here. 
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